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ABSTRACT 

A detailed analysis is presented of gravitational microlensing by intervening compact objects of the black hole shadows imaged 
by the Event Horizon Telescope (EHT). We show how the center, size, and shape of the shadow depend on the Einstein angle 
relative to the true/unlensed shadow size, and how the location of the lens affects the shift, size, and asymmetry of the black hole 
shadow due to microlensing. Assuming a supermassive black hole (SMBH) casts a circular-shaped true shadow, microlensing 
can create an asymmetry of up to approximately 8%, which is twice the asymmetry caused by the SMBH’s spin and its tilt 
relative to us. Furthermore, the size can be enhanced by ~50% of the true shadow. Currently, the terrestrial baselines of EHT 
lack the resolution to detect microlensing signatures in the shadows. However, future expansions of EHT, including space-based 
baselines at the Moon and Ly, could potentially enable the detection of microlensing events. For Sgr A*, an event rate of 0.0014 
per year makes the microlensing phenomena difficult to observe even with space-based baselines for the stellar population in the 
stellar bulge and stellar disk for lens masses ~ Mo. Nonetheless, the presence of a cluster of 20,000 stellar-mass black holes in 
the central parsec of the Milky Way, expected to arise from dynamical friction acting on infalling stellar clusters, significantly 
boosts the event rate. Hence, continuously monitoring the shadow of Sgr A* could offer novel insights into the compact object 


population surrounding the Galactic Center. 
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1 INTRODUCTION 


Black holes are natural laboratories for studying the most extreme en- 
vironments found in strong gravity. Event Horizon Telescope (EHT) 
exploits the Very Large Baseline Interferometry (VLBI) facilities 
across the globe, forming an effective baseline of earth size and re- 
sulting in a resolution of ~ 25 yas at 1.3 mm wavelength. One of the 
most astonishing discoveries of EHT has been the direct imaging of 
the shadows of the supermassive black holes (SMBH) at the center of 
the Messier 87 galaxy (M87* (Akiyama et al. 2019a,b,c,d,e,f)) and 
of the Milky Way (Sagittarius A* (Sgr A*) (Event Horizon Telescope 
Collaboration & Akiyama 2022a,b,c)). 

The theory that a black hole in the backlight of a background uni- 
form source at infinity would cast a shadow of its event horizon was 
first advanced by Bardeen (Bardeen 1974). Subsequently, (Luminet 
1979) and (Falcke et al. 2000; Falcke & Markoff 2013a) developed 
a more realistic shadow model by considering the emission from 
the accretion disc surrounding the black hole. The average angular 
radius of the shadow for a rotating black hole is widely understood to 
be (5+0.2)GM/(Ds c”), where Ds is the source’s distance from the 
earth and M is the SMBH’s mass (see eg. (de Vries 2000; Takahashi 
2004; Bozza et al. 2006; Johannsen & Psaltis 2010c; Psaltis et al. 
2015; Johannsen 2016)). The standard deviation about the average 
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radius arises from the asymmetric shape of the shadow, which can 
occur via changes in spin value and spin axis tilt (for a review of the 
precise dependence, see (Cunha & Herdeiro 2018; Gralla et al. 2019; 
Perlick & Tsupko 2022)). It is of the utmost importance to measure 
the size and shape of the shadow precisely because, assuming the 
Kerr metric describes the spacetime around a spinning black hole, 
it offers a unique probe for determining the spin and mass of the 
core SMBH. Furthermore, when the mass and spin of the SMBH 
are precisely determined through other observational methods, such 
as pulsars in proximity to the SMBH (Liu et al. 2014), stellar or- 
bits encircling the SMBH (Do et al. 2019), and the electromagnetic 
spectrum from accreting gas around the SMBH (Bambi 2017b), the 
observed characteristics of the shadow will offer a chance of identi- 
fying potential deviations from the Kerr hypothesis. This presents an 
opportunity to assess and test modified theories of general relativity. 


Another indirect technique for measuring the shadows involves 
photometric microlensing, where gravitational bending around mas- 
sive foreground objects leads to an apparent brightening in the back- 
ground source (Refsdal 1964; Irwin et al. 1989). This technique 
utilizes the microlensing phenomena induced by stars and compact 
objects within a galaxy that strongly lens a quasar, causing time- 
varying brightening (Kayser et al. 1986; Schneider & Weiss 1987; 
Kochanek 2004; Schmidt & Wambsganss 2010). It has been used to 
study accretion disks around the quasars (Grieger et al. 1988; Gould 
& Miralda-Escudé 1997; Agol & Krolik 1999; Fluke & Webster 
1999; Mortonson et al. 2005; Poindexter et al. 2008) and it has been 
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proposed as a tool for imaging the vicinity of the SMBHs in quasars. 
Although the quasar’s SMBH silhouette itself is too small for direct 
resolution, its accretion disk with a shadow creates a distinctive signal 
in the quasar’s light curve (Abolmasov & Shakura 2012; Mediavilla 
et al. 2015). This technique has also been extended to self-lensing in 
binary SMBH systems (D’ Orazio & Di Stefano 2018, 2020; Ingram 
et al. 2021). Furthermore, their shadows may self-lens due to strong 
gravitational effects and hence produce a dip in the light-curve (Dave- 
laar & Haiman 2022b,a), and a similar study has also been conducted 
in a stellar-mass black hole with binary companion star (Gott et al. 
2019). These scenarios of indirect shadow imaging hold consider- 
able promise for future observational endeavors (Falcke & Markoff 
2013b). 

Here, we present the formalism for the microlensing of the directly 
observable shadow. We apply the framework of microlensing on 
each point source of the boundary of the shadow due to a point 
lens in the foreground. Varying angular separations of the boundary 
sources from the lens lead to the apparent brightening and angular 
shift of each source on the boundary. This will be manifest as an 
overall distortion in the shadow. We will call the distorted shadow 
the microlensed shadow. In the case of relative motion between the 
source and the lens, the proposed distortion in the shadow will be a 
time-varying phenomenon. 

We present an analytic dependence of the microlensed shadow’s 
shape on the lens size and its variation with lens angular separation, 
assuming a circular shape of the true shadow. The primary signature 
of microlensing on the observed shadow is to enhance the size of the 
shadow and make the shadow asymmetric. The maximum enlarge- 
ment in the size of the shadow turns out to be 50% of the true size of 
the shadow. The highest asymmetry, however, can reach up to ~8% 
of the average radius of the microlensed shadow, which is twice the 
maximum asymmetry (4% (Psaltis et al. 2015)) arising due to incli- 
nation and the spin of the SMBH. As the non-zero spin of the SMBH 
may also cause asymmetry in the true shadow, a non-zero spin may 
contaminate the effect of microlensing. In this work, we restrict our 
analysis to the circular shape of the true shadow. We demonstrate the 
formalism of calculating microlensed shadow to the case of Sgr A* 
for which the abundances of stellar lenses in the Milky Way are well 
constrained observationally. We discuss the detectability assuming 
the achievable resolution from various baseline configurations such 
as Earth-size baselines (already achieved by current EHT observa- 
tions) and futuristic Earth-space baselines (Roelofs et al. 2019a; Fish 
et al. 2020; Mikheeva et al. 2020; Gurvits et al. 2022; Tiede et al. 
2022; Chael et al. 2023; Rudnitskiy et al. 2022). We compute the 
distribution of event rates throughout the microlensing event dura- 
tions by employing an analytic stellar model of the Milky Way. Our 
analysis indicates that the total event rate cannot exceed 107? per 
year. Nevertheless, we discuss scenarios where the event rate might 
experience enhancements. 

In a broader context, the ongoing microlensing event during the 
observation of the shadow can mimic the deviation in the size and 
shape of the shadow from the null hypothesis of the Kerr back- 
ground metric, which is also predicted by studies on searching beyond 
standard phenomena, such as testing the no-hair theorem of black 
holes (Johannsen & Psaltis 2010a,b), searching for super-spinning 
black holes (Bambi & Freese 2009; Bambi et al. 2019), testing the 
Kerr hypothesis (Bambi et al. 2012; Bambi 2017a) and probing light- 
scalar hair around the SMBH (Cunha et al. 2015, 2019). Therefore, 
before interpreting the deviation of a shadow from the null hypothesis 
in search of the beyond-standard phenomenon, one should rule out 
the standard possibility of a foreground compact object microlensing 
the shadow. Hence, the microlensing of a black hole shadow should 
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be treated as a standard astrophysical background phenomenon for 
the search of beyond standard physics using the shadow, which can 
be reduced by proper modeling of the microlensing event rate distri- 
bution over the event durations discussed in this paper. 

In sec. 2, we briefly review the microlensing of a point source due 
to a point lens. We then develop the formalism for microlensing of 
the shadow and give a generic expression for the distorted shadow 
caused by microlensing. In sec. 3, we characterize the microlensing 
signal of the shift in the shadow center and the radial variation for 
the microlensed shadow with position angle. The enlarged shadow 
will be described in sec. 3.3, and the induced asymmetry will be 
covered in sec. 3.4. Finally, in sec. 4, we discuss the detectability 
of the microlensed shadow of Sgr A* in an EHT extension with 
baselines at Moon and Lz positions. The anticipated duration of the 
events, optical depth, and the distribution of event rate over the event 
durations are calculated assuming a stellar population for the bulge 
and the disk of the Milky Way. The potential ways to enhance the 
event rate are also discussed. We then summarize our findings in 
sec. 5. 


2 MICROLENSING OF A BLACK HOLE SHADOW 


In this section, we will study the distortion of a black hole shadow 
due to the microlensing effect of a lens anywhere between Sgr A* and 
us. The shadow of Sgr A* is formed at the center of the Milky Way, 
where we have a huge stellar density in the central nuclear star cluster 
according to ref. (Gallego-Cano et al. 2018; Schédel et al. 2018; 
Baumgardt et al. 2018). For example, the total luminous stellar mass 
within 0.01 pc (1 arcsec) from Sgr A* is 180+30 Mo. Furthermore, 
the stellar bulge population (Valenti et al. 2016) of mass 2 x 10!° Mo 
within 1 kpc may also perturb the shadow. Hence, we can expect a 
significant perturbation of the shadow from the stellar components. 
Furthermore, the possibility of an orbiting IMBH (Deme et al. 2020) 
and stellar-mass black hole cluster in the central parsec of the Milky 
Way (Miralda-Escude & Gould 2000; Hailey et al. 2018) may also 
distort the shadow. The actual effect would depend upon the mass 
of the lenses as well. We will study the effect of lenses with mass 
ranging from planetary mass M ~ 1076 Mo and up to lenses of mass 
M ~ 104 Mo. 


2.1 Black hole shadow and lensing ring 


The shadow is one of the most distinguishing features of the direct 
image of an SMBH with a surrounding accretion disk. The photons 
from the accretion disk undergo strong gravitational lensing, which 
causes the shadow to form. Therefore, the “observed shadow” ex- 
tends up to the lensed position of the inner edge of the accretion 
disk (Luminet 1979; Falcke et al. 2000; Johannsen & Psaltis 2010c; 
Falcke & Markoff 2013a). 

Apart from the shadow, another feature of the black hole image is 
the lensed image of the optically thin accretion disk, which is yet to be 
observed by EHT or EHT-like VLBI facilities. The existence of the 
null geodesics surrounding black holes, whose order can be identified 
by the number of half revolutions of the photons around the black hole 
before escaping, causes the lensed images to develop. We will call all 
these lensed images “lensing rings” whereas some refer to all lensing 
rings with half rotations greater than two as photon rings. In principle, 
an infinite number of lensing rings exist, but all are formed around the 
critical curve (Bardeen’s shadow). Their thickness becomes sharper, 
and they form nearer to Bardeen’s shadow as the order of the lensing 
ring rises (Bozza 2010; Johnson et al. 2020; Gralla & Lupsasca 
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Figure 1. Ray diagram of gravitational lensing of a point source due to a point 
lens. The light rays coming from the source bend around the curved space 
around the lens (black-blob) and forms two apparent lensed images (Image 1 
and Image 2) of the source. 


2020). As a result, the size and shape of the lensing rings approach 
Bardeen’s shadow. They are thus asymptotically wholly governed by 
the characteristics of the black hole alone. Therefore, investigating the 
impact of the microlensing signal on Bardeen’s shadow is a natural 
choice that will not be affected by the accretion physics!. However, 
this paper’s formalism of microlensing of a shadow is equally valid 
for both the observed shadow and the lensing rings. We shall use the 
term boundary of the shadow to refer to any cases of the boundary 
of Bardeen’s shadow, the boundary of the observed shadow, and the 
boundary of the lensing rings without losing the generality of the 
effect. We will also assume a circular shape boundary of the true 
shadow. 

We will first review the microlensing of a point source due to a 
point lens. We then develop the theory of the effect of microlensing 
on the shadow. 


2.2 Brief review of microlensing 


Consider a point source in the sky at Ds distance away from us as 
shown in fig. 1. A point lens of mass M is also present between the 
source and us. Hence, its distance D; from us is always less than 
Ds. Due to the gravitational lensing effect, the light rays from the 
source would bend around the lens, giving the illusion of two separate 
point images of the source. To comprehend the geometry governing 
the production of the images, we imagine a lens plane that travels 
through the lens and is perpendicular to the line connecting the lens 
and ourselves. We can use the line between the source and us to 
project the source onto the lens plane. An angular separation vector 
B of the source with respect to the lens can then be defined such that 
Dı B represents the position vector of the projection of the source on 
the lens plane. The angular position vector of the two images can be 
ascertained by 


7 B 
= () 


6. = {1+ 11408 


l In principle, the shape of the lensing ring does not depend upon the ac- 
cretion physics, their observed shape may depend on it. This is because the 
accretion disk may have an inherent variability, and moving hot spot features, 
which can make the observed black hole image asymmetric (Akiyama et al. 
2019e; Zhu et al. 2019) and hence their lensing rings (Paugnat et al. 2022). 
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Figure 2. This plot depicts the relevant parameter space of mass M and radial 
distance r of the lens from Sgr A* towards Earth. Parameters in the colored 
region can potentially induce measurable distortion in the shadow expected 
to be ~ 0g, where the corresponding Einstein angle can be read from the 
color bar. Assuming 1% accuracy in the true radius Rshadow = R of the 
shadow, we obtain the lower cut-off in the parameter space by 8g = 0.01 R 
i.e. M ~ 6 x 1078 Mo(8200 pe/r — 1) for Sgr A*. This parameter region 
will be further refined later in the paper after quantifying the exact distortion 
due to microlensing. 


where 0g is the Einstein angle, and determines the angular scale of 
the separation between the two images and its dependence on the 
parameters of the source and the lens, scaling as (Gaudi 2012) 


2.85 mas EET (Be) (2) 


= IMo Ds \D; 

According to eq. 1, the angular position vector with respect to the 
lens is along B for one image and opposite to B for the other. As a 
result, the two images will always be aligned on the lens plane along 
the line connecting the lens and the projected source. There exists 
an imaginary ring of angular radius 6¢ centered at the lens called 
an Einstein ring such that one of the images would be outside of the 
ring, and the second image would be inside the ring. 

Another consequence of lensing is that the apparent brightness of 
the two images gets magnified as compared to the true brightness of 
the point source. The respective magnification in the brightness for 
the two images is given by, 

2 + 207, 
H = Z l (3) 


2662+402 7 


More details of the above discussion can be found in (Schneider et al. 
1992; Narayan & Bartelmann 1996; Wambsganss 1998). 

If the separation between the two images is small compared to a 
telescope’s resolution, then the two distinct images will not be ob- 
servable. However, the telescope would see a magnification-weighted 
average image called the centroid. The angular position of such a cen- 
troid relative to the lens will be (Hog et al. 1995; Walker 1995; Gould 
1996; Dominik & Sahu 2000; Nucita et al. 2017) 


ő = 0 la| +0 -lu-] _B es (4) 
c = i 
Ju+| + || B +262, 
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The centroid always lies on the line joining the lens to the projected 
source on the lens plane. The astrometric angular shift in the centroid 
relative to the true angular position of the source is quantified to be 


2 
Or, = 


Oc = ss 
2 2 
pr +20% 


be -B= (5) 
which is an angular vector defined such that Dz 6 c is the shift of 
the centroid on the lens plane from the true position of the projected 
source. The photometric magnification in the flux of the centroid 


would be, 


2 2 
pr +20, 


By B? + 407. l 


In the case of a relative proper motion of the source relative to the 
lens, we expect B to be time-varying; therefore, the angular position 
of the centroid would change with time. The phenomenon of variation 
of this apparent position of the source with time is called astrometric 
microlensing. On the other hand, a time variation in the magnification 
of the centroid is called photometric microlensing. 

Given that the shift in the centroid of a point source is ~ 80g, we 
can expect the distortion on a shadow of angular radius R would 
also be of the order of the Einstein angle. Therefore, assuming 1% 
accuracy in the measurement of the shadow, in fig. 2, we extract a 
target parameter space of the lens for which the Einstein angle and 
hence the distortion could be larger than 1%. 


He = lulta- (6) 


2.3 Theory of the microlensed shadow 


We will now summarize our formalism for describing the effect of 
microlensing on the boundary of a black hole shadow. We expect 
massive objects such as stars, stellar black holes, neutron stars, white 
dwarfs, brown dwarfs, dark matter, and even primordial black holes 
to surround a supermassive black hole at the center of a galaxy. 
Furthermore, these massive objects could also be floating between 
us and the SMBH. If these objects are close enough to the shadow 
along our line of sight during the observational period of the EHT- 
like telescope, they might serve as lenses for each point source of the 
boundary of the shadow. We expect that light rays originating from 
different points of the boundary of the shadow would get deflected 
due to gravitational lensing by intervening objects. 

To understand the effect of lensing on the boundary of the shadow, 
we assume that the true boundary of the shadow is circular. We also 
assume the boundary to be a collection of point sources. Although 
the influence of the lenses on each point source of the boundary 
will undoubtedly produce two images separated (~ 0g) from one 
another, it would be difficult for EHT-like VLBI facilities to resolve 
such effects. Hence we make an assumption that only the centroid of 
each point source would be observable by the telescope. We apply 
the formalism of microlensing described in sec. 2.2 to each point 
source that makes up the boundary of the shadow. 

In fig. 3, we depict a schematic diagram of the true boundary of the 
black hole shadow, represented by a red-solid circle. The foreground 
lens L is indicated by the black-filled blob on the lens plane. The 
true boundary is also projected on the same plane. The projected 
boundary of the shadow is shown as a orange-dashed curve on the 
lens plane. It is assumed that the lens is distance D; away from us. 
In contrast to the previous sec. 2.2, we now have numerous point 
sources as a part of the true boundary at the same distance Ds. 
However, the angular separation vector B of the point sources on the 
projected boundary from the lens would be different. Due to the small 
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Figure 3. This schematic cartoon illustrates the anticipated outcome from the 
microlensing of the shadow. The red-solid ring shows the true boundary of 
the shadow, which may be thought of as a collection of several point sources. 
The orange-dashed circle indicates the boundary projection onto the lens 
plane. As depicted in the cartoon, each point source on the boundary will be 
at a different angular separation from the lens L. Hence, it is expected that 
the shift in the centroid (given by eq. 5) will differ for various point sources 
across the boundary. The centroid of each point source on the boundary can 
then be located to determine the microlensed shadow. The expectation is that 
the microlensed shadow would be a closed contour of center, size, and shape 
different from the true boundary of the shadow. A cartoon of the microlensed 
shadow is shown with a red-dashed contour, which is distorted compared to 
the true shadow. 


angular scales involved, we can define a Cartesian coordinate on the 
lens plane with the center of the projected shadow as the origin. The 
X-axis is assumed to pass through the center of the projected shadow 
and parallel to the motion of the lens. Each point source is identified 
on the projected boundary by a position angle ø measured with the 
x-axis. The shift in the centroid due to the microlensing will not be 
the same for all the point sources because of their different angular 
separation from the lens. As a result, we can anticipate a distorted 
boundary of the shadow. The red-dashed curve shows the distorted 
shadow in fig. 3, which we will refer to as the microlensed shadow. 

We define the true boundary of the shadow to be a circular ring 
with an angular radius R to calculate the impact of microlensing on 
the shadow. The angular position vector of a point source on the 
boundary of the ring can be labeled by 


R= a 


R sin ġ D 
where ¢ € [0,27). Additionally, the angular position vector of the 
lens é is defined such that D; č is the position vector of the lens on 
the lens plane from the origin. As a result, the angular position vector 
of a point source on the true boundary relative to the lens would then 
be determined by 


R-¢, (8) 


and the magnitude of the angular separation between the lens and any 
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Figure 4. As an illustration of the phenomenon of microlensing of the black hole shadow, we have assumed the true shape of the boundary of the shadow of 
Sgr A* (Ds = 8.2 kpc) to be a perfect circle of radius R = 24.35 uas. The blue-dashed circle represents the true boundary of the black hole shadow. A lens is 
indicated by a black dot with an Einstein angle 0g = 7.8 uas. This Einstein angle corresponds to the mass of the lens to be M = 1 Mo and distance r = 0.5 
pe from Sgr A*. The microlensed images of an arbitrary source on the true boundary of the shadow trace the dotted contour which is the final microlensed 
shadow under centroid approximation (see sec. 2.3.1 for the validity). The color of the dots represents the photometric magnification of the images, where the 
magnification value can be read from the color bar. The shifted center of the apparent shadow is shown with a red marker. Various panels show the variation of 


the microlensed shadow as the lens passes in the foreground of the true shadow. 


point source on the boundary of the shadow is given by the following 
formula, 


Bp = Je + R2 - 2£R cos ¢. (9) 


Given the position of the point sources parameterized by ¢ on the 
boundary relative to the lens by eq. 8, we can obtain the apparent 
lensed location of the sources or centroids using eq. 4. The angular 
positions of these centroids expressed in terms of the position angle 
ġ are as follows, 


BA +367, = 
— = BU). 10 
B()? + 202, B(O) (10) 


The above expression of the centroid is an angular position vector 
from the lens location since eq. 8 provides the angular position of 
the centroid, which assumes the lens to be at the origin. But to obtain 
the microlensed shadow, the origin (denoted by O) has been chosen 
at the center of the true shadow. Hence, the locations of the centroids 
Ie with the true center as the origin are determined as 


6 (9) = 


fa. aD 
BCO)? + 20%, 
To obtain the expression of the contour of the microlensed shadow, 
we must first determine the new center of the apparent microlensed 
shadow labeled by C. We define the shifted geometrical center of the 
microlensed shadow as 


Olc(d) = 6c+E=R+ 


OÏ e(n) + OF (0) 
5 : 

We will call this the lensed center. The lensed center of the mi- 
crolensed shadow will always align with the lens and the center of 
the true shadow. A detailed dependence of the lensed center on the 
lens and true shadow characteristics will be discussed in sec. 3. 

Now, to quantify the shape of the contour of the centroids with 
regard to the boundary of the microlensed shadow, we introduce an 
angular radial vector R L(@). It starts at the lensed-center and ends at 
the centroid of a point source on the true boundary. The expression 
of R L(¢@) originated from the shifted center C is given by, 


OC = (12) 


; : po f 62 z 
Rr =0l.-0C=R-0C+ E P(g). (13) 
E 


B(g)? +20 


It should be noticed that the position angle of R L() will not always 
be the same as the position angle (¢) of R. This suggests that the 
microlensing will also cause a twist along the boundary. The angle, 
say L, formed by R L(@) with the positive x-axis can be found by 
using the components of R L(@) along the x-axis and the y-axis as 
indicated by Ry, and Rz ,y, respectively. Hence, the twisted position 
angle is obtained by, 


R 
iy ss: See Poi (14) 


Rix 


Another consequence of microlensing on the shadow would be 
non-uniform brightness magnification for the various portions of the 
boundary. The magnification in the flux (see eq. 6) varies for the 
point sources across the boundary as a result of the same fact that 
the angular separation £ differs for different sites on the boundary. 
Therefore, the non-uniform brightness of the boundary would be a 
sign of photometric microlensing of the shadow. A drawback of using 
eq. 6 is in the case of the lens lying on the ring. The point source 
perfectly aligning with the lens will have the magnification to be 
theoretically infinite, but practically will be limited by the total energy 
emitted by the point source per unit time and per unit area around 
the point source, and hence remains finite. Although the expression 
of magnification under the point source approximation blows up for 
the perfect alignment, the magnification-weighted average angular 
position of the image (centroid) given by eq. 10, which is used to 
calculate the distorted shape of the shadow, will be zero for the point 
source perfectly aligned with the lens and remains finite for all other 
values of the lens-source separation. 

In fig. 4, we demonstrate the effect of microlensing on the true 
boundary of the shadow shown by the dashed-blue circle having a 
center at O. To illustrate the effect, we consider Sgr A* located at a 
distance Ds = 8.2 kpc. We assume the true black hole shadow has 
a circular shape with a radius of R = 24.35 yas. The figure exhibits 
how the microlensed shadow changes for various angular positions 
of the lens of mass 1 Mo and at a distance of r = 0.5 pe from 
Sgr A* in the foreground of the shadow. Here, we display the true 
boundary of the shadow, the true center, and the lens, together with the 
microlensed shadow and the lensed-center. These panels reveal a few 
general characteristics of the microlensing of the shadow. First, the 
shadow shifts away from the lens, which is shown by the lensed-center 
(see the red marker). However, when the lens and the true shadow 
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are perfectly aligned, the lensed center also coincides with the true 
center. Second, the shadow’s shape is warped, and its size is enlarged. 
The shadow takes on the shape of a circle when the lens is perfectly 
aligned with the true shadow. The third effect is the asymmetry in the 
brightness of the boundary. This is due to the uneven magnifications 
for various portions of the boundary of the shadow. In the case of the 
lens and the true shadow being perfectly aligned, the only remaining 
characteristics of microlensing are the shadow’s enlarged size and an 
overall magnification in the brightness of its boundaries. However, 
the brightness asymmetry and the asymmetric shape vanish when 
the lens and the true shadow are perfectly aligned. 

The figure shows that the boundary’s asymmetric brightness is 
important when the lens is close to the true shadow. In contrast to 
the asymmetry in brightness, the asymmetric shape of the shadow 
survives with a wider angular separation of the lens from the true 
shadow. This behavior can be understood qualitatively using a point 
lens’s microlensing of a point source. It is well-known that the photo- 
metric magnification signal fades away much more quickly than the 
astrometric shift signal at large angular separations of the lens. This 
is because the astrometric shift (given by eq. 5) falls off as ~ 1/2 at 
large 8, whereas the photometric magnification (provided by eq. 6) 
falls off as ~ 1/64. That is why, even if the lens were far from 
the shadow, the strength of the astrometric distortion is substantially 
stronger than the photometric distortion. 

The photometric asymmetry caused by microlensing could be im- 
portant for the event when the lens passes very close to the boundary. 
It will change the morphology of the boundary, and microlensing 
being a time-varying phenomenon, it will add another time-scale 
(Einstein crossing time tg = @¢/p, where p is the angular velocity 
of the lens relative to the source) into the observability of the bound- 
ary of the shadow apart from the innermost stable circular orbital 
time scale. However, if the integration time (time taken by EHT-like 
telescope to collect enough photons to calculate the visibility) is 
larger than tg, the photometric effect will be washed out. Given that 
the astrometric distortion will have its influence from a much larger 
angular separation of the lens from the shadow as compared to the 
photometric effect, we will be focusing on the impact of astrometric 
microlensing on the shadow. 


2.3.1 Centroid approximation and its validity 


It should be noted that fig. 4 shows the microlensed shadow under 
the centroid approximation. In principle, each point source on the 
boundary of the true shadow will have two lensed images due to 
the point lens as discussed in sec. 2.2. Hence, there will exist two 
lensed images of the shadow at any instant. In fig. 5, we illustrate 
two lensed images of the shadow, color-coded to represent their cor- 
responding photometric magnifications, along with the microlensed 
shadow under the centroid approximation outlined by the red con- 
tour. The lens and source parameters are still same i.e. Ds = 8.2 kpc, 
R = 24.35 pas, M = 1 Mo, andr = 0.5 pc. 

In the left panel of the figure, we consider a scenario where the 
lens (depicted as a black dot) is sufficiently distant from the true 
shadow (indicated by the blue dashed circle) such that the angular 
separation (£) between all points on the true boundary is larger than 
the Einstein angle (8g), i.e., 8 > 0g. Under these conditions, the 
primary shadow to the right of the lens exhibits an approximately 
unit photometric magnification throughout its extent. Conversely, 
the secondary image, positioned to the left of the lens, experiences 
significant demagnification. This distinction arises from the fact that 
the photometric magnifications of the primary and secondary images, 
as given by eq. 3, exhibit variations described by u+ ~ 1+ OF / pt 
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Figure 5. The figure displays two lensed images of the shadow of Sgr A*, 
with contours color-coded to represent the corresponding photometric magni- 
fication. The red contour represents the centroid approximation. The lens and 
source parameters are consistent with those in fig.4. In the left panel, the an- 
gular separation of the lens is large, resulting in significant de-magnification 
of the secondary lensed image, leading to the primary lensed image predom- 
inantly contributing to the measured flux. The primary lensed image closely 
matches the centroid-approximated contour in this scenario. In the right panel, 
the lens is close enough for a portion of the shadow boundary to lie within 
0.50g around the lens. Under such conditions, the secondary image becomes 
observable if the angular resolution is finer than 6g. 


and u- ~ OF f Bt respectively when £ is much larger than 0g. When 
observing such a lensing scenario of the shadow with a telescope 
like EHT, the flux contribution from the secondary image becomes 
negligible. Consequently, the primary lensed shadow, to which the 
centroid shadow is approximated, emerges as the observable feature 
of interest. 

On the other hand, it can be calculated from eq. 3, when the lens 
is close enough to the shadow such that there is a portion of the 
true shadow boundary for which £ < 0.50g, viz. shown in the right 
panel, then the secondary image of that portion can be magnified and 
roughly would be separated from the primary shadow with ~ 0g. 
Hence, the secondary shadow can, in principle, be observable along 
with the primary shadow, except the telescope resolution per epoch 
is less than the Einstein angle. 

In summary, the centroid approximation is valid as long as the Ein- 
stein angle is smaller than the resolution per epoch of the telescope. 
However, when the Einstein angle is larger than the resolution per 
epoch, the observability of an additional secondary shadow, along 
with the distorted primary shadow, necessitates the lens to be po- 
sitioned sufficiently close to the true shadow such that 6 < 0.50g 
for some point sources along the true shadow boundary to remain 
magnified (or at least not demagnified). Nevertheless, even with sub- 
Einstein-angle resolution, if 8 >> 8g, the highly demagnified sec- 
ondary shadow can be disregarded, and the centroid image can be 
accurately approximated to the primary shadow. Therefore, to de- 
termine the maximum allowable £ for which the distortion is mea- 
surable, the centroid approximation suffices. Consequently, for the 
purpose of studying the detectability of this phenomenon, we will 
continue to employ the centroid approximation throughout the paper. 

In the next section, we will characterize the astrometric microlens- 
ing signal imprinted on the boundary of the shadow. 


3 CHARACTERIZATION OF THE MICROLENSING 
SIGNAL 


We first define the signals of microlensing of the shadow and then 
we study their dependence on the lens and the true shadow. 


3.1 The lensed center 


Equation 12 provides the location of the lensed center OC. If the 
angular separation of the lens from the true center is € and the 
Einstein angle is 0g, then the expression of the location of the lensed 
center will be 
2 2 2 2 
We a: 07, (205, -R + &*) a (15) 
[207 + (R - €)7] [20]; + (R + £7] 


The numerator of the above expression suggests that in the case of 
the radius of the true shadow smaller than or equal to VŽ times the 
Einstein angle (i.e., R < V20g ), the shift in the center will be along 
the negative x-axis for all £ > 0. However, in the case of R > V20Ẹ, 
the shift in the center will be along the negative x-axis for the angular 


separation greater than equal to ,/R2 — 262; otherwise, it would be 
along the positive x-axis. The lensed-center will coincide with the 
true center for a perfect alignment (€ = 0) of the lens with the true 
shadow. This is because, for perfect alignment, each point source on 
the boundary of the shadow shifts symmetrically; hence, the average 
position of the point sources at the boundaries would be unaltered. 

In fig. 6, we plot the magnitude of OC, denoted by OC in the 
figure, with £. A negative value of OC indicates that the lensed 
center has moved away from the lens due to microlensing. On the 
other hand, a positive value of OC shows the lensed center is on the 
same side as the lens from the true center. A vertical black-dashed 
line marks the true boundary of the shadow, which is assumed to 
be of Sgr A*, at R = 24.35 yas. This figure suggests two important 
lessons about the location of the center of the apparent shadow. First, 
if the radius of the true shadow is smaller than V29F, the lensed 
center will always be on the opposite side of the lens from the true 
center. In contrast, there will be a sign change in OC if the radius of 
the true shadow is greater than V20g. The sign change in OC occurs 
near the boundary of the true shadow, which is evident from eq. 15. 
Intuitively, this behavior becomes straightforward when we recall 
that the astrometric shift of any point source at the true boundary 
of the shadow is always away from the lens. This shift increases 
as 8 increases, reaching its peak at 8 = V28g, and then gradually 
decreases with further increases in 8. Consequently, a lens located 
outside the true shadow will consistently push all the point sources 
on the boundary away from the lens, causing an overall shift of the 
shadow away from the lens. This holds true regardless of the size of 
the shadow relative to the Einstein angle. However, for a lens within 
the true shadow, the shift can be towards the lens when the closer side 
of the true boundary shifts more than the far side, which is possible 
only when R < V28g. Hence, for lenses with a tiny Einstein angle 
as compared to R, the lensed center will be on the opposite side of 
the lens from the true center if the lens is outside of the true shadow, 
whereas it will be on the same side of the lens if the lens is within 
the true shadow. Second, there are angular positions of the lens for 
which there will be an extreme shift in the center. It can be seen from 
the figure there will be two é values for which shift attains extremum 
values if Einstein angle is smaller than R (specifically R > V20R). 
One extremum is attained when the angular position of the lens is 
outside of the true shadow, and the other occurs when the lens’s 
angular position is inside the true shadow. 

We also draw the same curve for various Einstein angle values 
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relative to the shadow’s true radius. The vertical black-dotted line 
represents the true boundary of the shadow at € = R. In the case 
of Og > R/V2, typically, the shift grows with increasing angular 
separation. However, it attains a maximum value when the angular 
separation is around 26, and later as the angular separation is 
larger than this value, the shift goes down. The shift scales as 02, [£ 


for large angular separation (£ >> 0g). Inthecaseof 8g < R/V2, the 
shift has a non-trivial dependence near the boundary of the shadow. 


It changes sign near € = ,/R2 — 202, evident from the expression 
given by eq.15 as well. 

It should be noted that the absolute shift in the center of the shadow 
is not an observable signal in a single epoch of the image. Because if 
the shadow is microlensed, then the true center would not be observ- 
able. Therefore, a single frame of the shadow captured at any epoch 
would be unable to measure the shift in the center due to microlens- 
ing. However, one can observe the change in the shifted center if at 
least two frames of the shadow are taken at two different epochs. 
Therefore, the change in the position of the lensed-center between 
the two frames could, in principle, be an observable lensing signal. In 
practice, accurate astrometry on the shadow center is crucial to mea- 
sure the shift in the center, which necessitates the use of a standard 
source in the field of view to act as a stable reference point. However, 
finding a well-calibrated source with precisely known astrometry to- 
wards the galactic center is extremely difficult due to highly crowded 
stellar density and high extinction, for eg. (Lindegren et al. 2021). 
Furthermore, in the case of the duration of the microlensing event 
being much longer than the observation time, even any change in the 
lensed center would be unobservable. 


3.2 Variation of shadow radius with position angle 


Another feature of microlensing of the shadow is the variation of 
the radius with the position angle. A lensing signal is defined as 
the variation of the lensed radius Rz, with respect to the observed 
position angle @, at any instant. The exact dependence of Rz (@) and 
oL(@), given by eq. 13 and eq. 14 respectively, on the Einstein angle 
and the angular separation is not trivial to get an overall physical 
understanding. Therefore, we break our explanation into two cases 
(1) a much larger Einstein angle compared to R and (2) a much 
smaller Einstein angle compared to R. 

Firstly, we take the case when Einstein angle is much larger than 
the true angular radius of the shadow, i.e., 9g >> R. The radial 


profile of microlensed shadow can then be approximated as 
RE z% Ro = Ao cos2¢ (16) 


where the mean radius Rp and the amplitude Ao of the oscillation 
are given by 


3 RR RE 
R-RE -RÉ se 26x, 
Ryo = {287 4a Dee E SOE (17) 
R ;é > OB; 
2 
RE ;E < OE, 
Ag = g (18) 
Ra 5€ > OE; 


and the observed position angle is approximated by 


Bava i aay. exe 
a Page nC) + ggz Snl ġ) é< 0E, ie 


Or, 
+ z3 sin(24) sé > OE. 
In fig. 7, we show how the exact variation of the angular radius with 
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Figure 6. This figure shows how the angular position of the lensed-center OC, relative to the center of the true shadow, varies with the angular separation of the 
lens with the true center. The center is supposed to be moved in the direction of the lens when OC has a positive value. However, a negative value of OC shows 
the lensed-center shifts away from the lens. The radius of the true shadow, which we have taken to be of SgrA* (R = 24.35 uas) for demonstrative purposes, 
is shown by the vertical black-dashed line. In the left panel, we plot OC for the cases when the radius of the true shadow is greater than V20g. Every time 


the value of £ crosses the value of 4/ R? — 202, the sign of OC changes. There are two extremum shifts in the center. One occurs when the lens is outside the 


shadow, and the associated extremum shift is away from the lens. The second occurs when the lens lies within the true shadow, and the corresponding extremum 
shift is towards the lens. For a small Einstein angle compared to the true size of the shadow, the shift in the center due to microlensing would be substantial only 
if the lens is close to the boundary of the true shadow. In the right panel, we present the cases for which the radius of the true shadow is less than or equal to 
V20Ẹ. It is important to note that, in this case, the shift will always be away from the lens, and now, the extreme shift occurs only when the lens is outside of 


the true shadow. 


the position angle changes as the lens-shadow separation € changes. 
The black-dashed horizontal line represents the true radius R. The 
oscillation in the radius about a mean radius, as the position angle 
changes from 0 to 27, is the primary characteristic of the subplots 
of the figure. We also show the dependence of the signal on the 
Einstein angle in the same figure. For a fixed £, it can be seen that 
the curve of Rz with ¢ with a larger value of 0g is always above 
the curve with a smaller value of 0g. It implies that the average size 
of the shadow would be larger for the larger value of the Einstein 
angle for a given lens-shadow separation. The above approximate 
expression of Ro captures this dependence explicitly. Furthermore, 
it is also evident from the dependence of Ro that the average radius 
saturates to a value of 3R/2 — R3/ (402) when the lens coincides 
with the true shadow shown in the subplot. Hence, microlensing 
cannot arbitrarily increase the size of the shadow with arbitrarily 
large 0g. In fact, if 0g >> R, then the maximum size attained by the 
microlensed shadow is 3R/2 which is independent of 0g. Therefore, 
microlensing can only magnify the shadow up to 50% of the true size 
of the shadow. 


In contrast with the size of the shadow, there are oscillations in the 
radial profile, which are a signature of apparent asymmetry in the 
microlensed shadow. The oscillation amplitude increases as the lens 
moves away from the shadow, starting from £ = 0 yas (e.g., see the 
amplitude of the red curves in the subplots). It is also worth noticing 
that in the far region around the shadow, the oscillation amplitude 
also goes down as the lens asymptotically moves away from the lens. 
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This asymptotic behavior is well captured by the Ag expression given 
by eq. 18. It suggests that there will be an angular separation of the 
lens from the shadow for which the amplitude of the oscillation in the 
radius will be maximum and hence the asymmetry will be maximum. 
We numerically calculate this distance to be ~ 26g in the limit of 
large Einstein angle”. The maximum amplitude Ao,max turns out to 
be ~ R/8 when £ ~ V26p. 

Secondly, we consider the case when the Einstein angle is much 
smaller than the true radius of the shadow, i.e., 0g << R. The 
microlensed shadow profile can again be approximated by 


Rg =œ% ro- agcos2ġ (20) 


where the mean radius rg and the amplitude ag of the oscillations 
are given by 


Rete 2% Seek 

ro = RIR ESR (21) 
R 5E>R; 
OE? 
- é< R, 

a = y (22) 
RÆ ;R <Æ; 


E 


2 Notice that the shift of a centroid of a point source is maximum when the 
angular separation is V20E 


Microlensing BH Shadows 9 


& = -1000 pas & = -750 pas = = -500 pas & = -250 pas = = 0 pas = = 250 pas & = 500 pas = = 750 pas = = 1000 pas 
ININ 
1 ANIN 1 1 
0 In 2n 0 In 2n 0 in 2n 0 in 2n 0 in 2n 
9 % 9 $ 9 
ININ 
1 1 1 1 1 | 
0 in 2n 0 in 2n 0 in 2n 0 in 2n 0 1n 2n 
$L $L $L $L $L 
L 1 L 1 L L 1 1 
0 in 2n 0 in 2n 0 in 2n 0 in 2n 0 in 2n 0 1n 2n 0 in 2n 0 in 2n 0 1n 2n 
$ $ $ 9 $ $ $ $ 
— 6¢= 100 pas —— 6 = 200 pas — 06¢ = 300 pas — 6¢ = 400 pas —— e= 500 pas 


Figure 7. We describe the profile of the angular radius of the microlensed shadow with position angle. Here we show the angular radius profile when the Einstein 
angle is larger than the angular radius of the shadow. The horizontal dashed-line shows the true angular radius of Sgr A*. Each column corresponds to different 
angular separations of the lens with the center of the true shadow. The angular separation reduces from left to right, coincides with the center of the shadow at 
the central panel, and increases again. In the first row, we plot the observed angular radius with the position angle of the true shadow. The second row shows 
the observed profile of the angular radius of the microlensed shadow with the observed position angle. However, the third row is drawn to show the difference 
between the apparent position angle and the true position angle of a point of the shadow and how it changes with the true position angle. The typical profile is 
well explained by eq. 16 and eq. 19. 
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Figure 8. This figure shows the variation of the angular radius profile of the microlensed shadow similar to fig. 7 but when the Einstein angle is smaller than 
the angular radius of the true shadow. Here we increase the angular separation across the columns from left to right. The basic feature of the profile can be 
encapsulated in eq. 20 and eq. 23 (for more details see the text). 


and the observed position angle is approximated to be 


ot 
ot 


L 7 


a sin(ġ) + 


sin(2¢) 


4 


2 92 
TE sin(2¢) 


E< R, 
R< È. 


values of €. Again the black-dashed horizontal line represents the 
true radius R. The mean radius rg is smaller for a smaller value of 
Einstein angle, and asymptotically it coincides with the true radius 
when 6¢ — 0. However, the amplitude of the oscillation ag increases 
as the lens moves away from the center of the shadow. Furthermore, 


(23) 


In fig. 8, we plot an exact dependence of the Rz on ġz for various 
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Figure 9. This figure depicts how the average angular radius changes as the 
angular position of the lens changes. The averaged angular radius approaches 
a maximum value as the angular separation of the lens is $ Oz. It also shows 
that microlensing can not arbitrarily magnify the shadow. The maximum 
enlargement of the shadow that can occur through microlensing is 1.5 times 
the size of the true shadow when the Einstein angle is much larger than the 
true shadow. 


for very large £, it again decreases. These behaviors of rg and ag are 
well described by eq. 21 and eq. 22 respectively. 


3.3 Magnification in size 


We define the radius of the microlensed shadow as the average of 
Rz (¢ġ) over . Hence, the average radius is given by 


1 2m 
(Ri) = =f Rola. (24) 
TJO 
In the limiting cases of 0g >> R and 0g << R, we have obtained 
an approximate expression of Rz in the previous sec. 3.2. Hence, the 
average radius in various limiting cases can be summarized as 


Ro ;0E >>R, 
(RL) = (25) 
ro 306 << R; 


where Rọ and rọ are given by eq. 17 and eq. 21 respectively. 

In fig. 9, we show the change in the average size of the microlensed 
shadow with respect to the true size of the shadow as a function of 
the angular separation of the lens é. It can be seen from the plot 
that the size of the microlensed shadow gets magnified as the lens 
comes closer to the shadow, but there is a maximum magnification 
in size. There is a plateau region as the lens-shadow separation is 
less than around V20F. In the limit of Og >> R and £ << Og, the 
leading order behavior of Rz — R would be Ro ~ R/2. Hence, the 
maximum change in the size of the shadow will be half of the true 
size of the shadow. However, in the limit of R >> 0g and £ << R, 
the plateau goes down as the Einstein angle reduces, which can also 
now be understood with the leading order behavior of the Rz, — R as 
ro~ 02, /R. The dependence of the change in the size of the shadow 
on the explicit parameters of the lens such as mass and the distance 
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of the lens from the center of the Milky Way can also be seen from 
fig. A2 in the appendix. 


3.4 Asymmetry in shape 


The oscillation in the radius versus position angle curve is a signature 
of the asymmetry in the shadow, as shown in fig. 7 and fig. 8. Unlike 
the shift in the center and the magnification in the size of the shadow, 
asymmetry is an absolute observable effect of the microlensing of 
shadow. Hence, it can be measured even in a single epoch frame of 
the shadow. The asymmetry of the lensed shadow is defined as the 
variance in the observed radius of the shadow and given by, 


1 27 
A = i af RL- Rd 28) 
z Jo 
The asymmetry in the various limits can be summarized by 
Ao. 
= ÂE >>R, 
A=% Phe one. (27) 
vy? E << K; 


where the oscillation amplitudes Ap and ao are given by eq. 18 and 
eq. 22 respectively. 

In fig. 10, we show the exact dependence of the asymmetry on 
the angular separation € and the Einstein angle 0g. The curves with 
different colors signify the various values of the Einstein angle. There 
are two distinct behaviors of the curve depending upon whether the 
true radius of the shadow is less than V28 E or greater than V20r. 

In the case of R < V20g, shown in the right panel of fig. 10, 
the asymmetry increases with an increase in €, reaching a maximum 
value at £ = V20g, and then goes down with a further increase 
in the angular separation é. This is because, for the large Einstein 
angle, the amplitude Ag is given by eq. 18 which for small £ varies 
as ~ é*/ 0%. However, when the lens is far from the shadow (£ > 
0p ), Ag ~ 03 JE? and decreases as we increase the lens-shadow 
separation. The approximate separation for which there will be the 
maximum asymmetry is obtained numerically to be ~ V20g and 
the corresponding maximum asymmetry ever achieved due to the 
microlensing effect is found to be 


R 

Amax aya" (28) 
which is independent of the Einstein angle itself. This explains why 
the maximum asymmetry in fig. 10 is the same for all 9g >> R/V2. 
On the other hand, two maxima occur in the case of R > V20E 
which is shown by the left panel of the figure. The first maxima 
occurs when the lens is within the true shadow, and another occurs 
outside. In this limit of a small Einstein angle compared to the true 
size of the shadow, the asymmetry signal would be significant when 
the lens is near the boundary of the shadow, also shown in the left 

panel of the figure. 


4 DETECTABILITY OF THE MICROLENSED SHADOW IN 
AN EHT-LIKE TELESCOPE 


We shall now estimate the telescope specifications needed to detect 
the microlensing of a black hole shadow. This requires quantifying, 
first, the uncertainty in measuring the shape of the shadow with a 
given telescope specification and second, how often a lens will be 
aligned with the shadow to induce a detectable signal. To demonstrate 
the procedure, we shall again concentrate on the microlensing effect 
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Figure 10. This figure depicts the variation of the asymmetry in the microlensed shadow with the lens angular separation, which is calculated numerically 
using the exact expressions given by eq.13 and eq.26. Different curves are drawn for different Einstein angles represented with different colors. We separate the 
asymmetry signal in the above two panels depending upon the relative size of the lens relative to the black hole shadow. The dashed-vertical line is the angular 
radius of the true shadow, which is again taken to be of Sgr A*. Left panel shows the case of Einstein angle less than R/V2, which approximately implies a 
smaller Einstein angle as compared to the angular radius of the true shadow. There are two maxima in the asymmetry, the first maxima occurs when the lens is 
outside of the boundary of the true shadow, and the second occurs when the lens is within the shadow’s boundary. For a sufficiently small Einstein angle, the 
asymmetry signal is mostly strong when the lens is near the true shadow’s boundary. Right panel is the asymmetry signal for a larger Einstein angle as compared 
to the size of the shadow (Og > R/V2). In this case, the maximum asymmetry is obtained for the angular lens position outside the true shadow. Asymptotically 
for a very large Einstein angle, the maximum asymmetry approaches the order of the radius of the true shadow. 


on Sgr A* for which the abundance of stellar lenses in the Milky Way 
are well constrained observationally. 

Apart from the shift in the center, the enlargement in the aver- 
age size, and the asymmetry in the microlensed shadow, we saw in 
the previous sec. 3 that the asymmetry is an absolute observable 
(assuming the circular form of the true shadow). Additionally, the 
average size of the shadow of Sgr A* is tightly constrained with 
precise measurements of its distance and mass by the stellar orbits 
around Sgr A* (Boehle et al. 2016); therefore, if the measured size 
surpasses the expected size from the stellar orbits, it indicates an 
observable effect of microlensing. The approximate expressions for 
the asymmetry and the average size up to the leading order in é are 
derived under the assumption that the lens is far from the boundary 
of the shadow compared to the Einstein angle. They are given by 


02 
A z e (29) 
and 
9 OF. 
(RL) =x a1 305) (30) 


Considering the uncertainty in measuring the shadow radius as Ca, 
which represents the standard deviation of the radial measurement, 
we determine the maximum threshold angular separation éh of the 
lens from the center of mass of Sgr A* such that A = oq and 


(RL) — R = Ga, as follows: 


19-3 \ 1/2 
ha = 26.60g (a) > (31) 
a 
for asymmetry and 
10-3 1/4 
hR = 07e (2 al (32) 
a 


for size magnification. These thresholds define the angles within 
which a lens must lie to induce a detectable microlensing event cor- 
responding to each observable. As the asymmetry observable dies 
out more slowly (67, /é2) compared to the size magnification (OF. /é4) 
with increasing lens angular separation &, the asymmetry observable 
will give us a maximum event rate due to the larger threshold angle 
compared to the size magnification. There will be a subset of the 
maximum number of events, which will show a detectable magni- 
fication in size. We will see the event rate scales linearly with éh, 
and the threshold separation for the two observables are related by 


Eth,A © 38 (10-3/(oa/R)) &th,R,,- Hence, we utilize the asym- 
metry observable as a proxy to determine the total microlensing 
event rate. The subset of events inducing detectable size magnifica- 
tion can then be obtained by dividing the total event rate with the 
factor 38 (10-3/(oa/R)) 

According to sec. 3.4, the asymmetry in the shadow due to mi- 


crolensing results in variance (02) in the radial measurements of the 
shadow. Regardless of how large the Einstein angle is, the max- 
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Figure 11. This figure shows the parameter space of the lens to be probed by various possible radio baselines in space by observing the shadow of Sgr A*. The 
color represents the asymmetry induced in Sgr A* due to microlensing by a lens, in our line-of-sight, of a given mass and distance from the center of the Milky 
Way. Various panels show how the asymmetry and hence the region of the parameter space changes as the angular separation é of the lens changes with distance 
from the shadow. The region above solid lines and dashed lines contours are the required characteristics of the lens for showing asymmetry signals in various 
baselines at Moon and Ly positions, respectively. The black and green contours correspond to 1.3 mm and 0.5 mm operating radio wavelengths, respectively. 


imum standard deviation in the radius of Sgr A* can be up to 
Amax = 2.152 uas due to microlensing. Therefore, the radius of 
the shadow cast by Sgr A* must be measured with a precision of at 
least 2.152 uas in order to have any microlensing signal. The current 
Event Horizon Telescope has measured the diameter of the shadow 
to be 2R = 48.7+7.0 uas (Event Horizon Telescope Collaboration & 
Akiyama 2022a). Hence, there cannot be any microlensing signature 
in the currently available image of Sgr A* by EHT. However, given 
the future potential of EHT and the planned space-based EHT-like 
VLBI facilities (Roelofs et al. 2019a; Fish et al. 2020; Mikheeva et al. 
2020; Gurvits et al. 2022; Tiede et al. 2022; Chael et al. 2023; Rud- 
nitskiy et al. 2022), we anticipate achieving direct imaging of Sgr A* 
with a sub-microarcsec level of resolution. As a result, it should be 
possible to observe the microlensing phenomena that we suggest via 
direct imaging of the SMBH shadow. 

We can estimate the projected precision in the radial measurements 
that could be made at EHT-like VLBI facilities in the near future. 
These estimates will allow us to find the range of mass and distance 
of the lens, which may result in an observable microlensing effect on 
the black hole shadow of Sgr A*. 

Let us assume that the largest baseline that can be used is D and 
that the radio observation’s wavelength is 4. This will enable the 
maximum resolution of the black hole shadow to be res ~ A/D per 
epoch. We consider an annular region with an angular radius of R 
and an angular thickness of Ores. Hence, the number of data points N 
on the boundary of the shadow is given by, 


2nROres _ 2R 


N = 5 = : 
TO res Ores 


(33) 


This implies that we will have N number of radial measurements of 
the shadow corresponding to the N directions ¢ in the interval of 
ġ € [0,27) with step size zm. Each radial measurement will have 
an accuracy of Ores. Therefore, assuming the true boundary of the 


shadow to be circular, the standard deviation of the angular radii 
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about the true radius R would be 


Ores 
Oa = — 


VN 


Q 


A 106 km\? 24. 
0 =) 35pas (34) 


0134 
0:013% pras (a D R 


The percent accuracy achieved will then be given by 


ferror = 


2 


3/2 
A 106 km 24.35 pas 
0.055% maa F ) . (85) 


We examine the detectability of the microlensing phenomenon 
in the black hole shadow of Sgr A* through three benchmark base- 
line configurations. First is Earth configuration, where the maximum 
baseline length corresponds to the size of Earth, akin to the configu- 
ration employed by the EHT. Second, Earth-Moon configuration con- 
siders a setup where a radio telescope is positioned on Earth, while 
another is situated on the Moon, creating a maximum baseline equiv- 
alent to the Earth-Moon distance. The third configuration Earth-L 
consists of a radio telescope located at L3 (the second Lagrange point 
of the Earth-Sun system) in conjunction with one on Earth, forming a 
maximum baseline separated by Earth-L>. In each configuration, we 
consider two operational radio wavelengths, namely 1.3 and 0.5 mm. 
For the Earth-Moon and the Earth-L2, we present the uncertainty in 
the radial measurements in each epoch because of the unknown ca- 
dence and imaging time (tim), where the imaging time is the required 
time to make one black hole image, which would depend upon how 
these arrays are implemented for these space-based baselines. How- 
ever, for the Earth-based arrays currently available in EHT operating 
at 1.3 mm, we know the observation scheme. According to the EHT 
collaboration (Event Horizon Telescope Collaboration & Akiyama 
2022b), there are two 8-minute scans of Sgr A* per hour, followed 
by ~3-minute scans, with an 8-minute gap. This observation scheme 


A= 1.3 mm 
D (km) Ores N # of Ta Jo 
(uas) epochs (uas) error 
10,700 (Earth) 25.06 1.9 1 17.98 73.8 
300,000 (Earth-Moon) 0.89 54.5 1 0.121 0.50 
1,500,000 (Earth-L2) 0.18 272.4 1 0.011 0.04 
A = 0.5 mm 
10,700 (Earth) 9.64 5.1 1 4.288 17.61 
300,000 (Earth-Moon) 0.34 141.7 1 0.029 0.12 
1,500,000 (Earth-L2) 0.07 708.3 1 0.003 0.01 


Table 1. This table shows the uncertainty in the measurement of the radial 
measurements at various baseline configurations. It is used to make predic- 
tions for the observability of the asymmetry in the microlensed shadow of 
Sgr A* due to the stellar granularity of the Milky Way. 


has been followed for about a day to obtain one image, i.e. the imag- 
ing time for EHT is currently about a day. We present the accuracy 
achieved by these configurations of the possible baselines in tab. 1. 

Using these benchmark resolutions of various space-based base- 
lines, we derive the in-principle range of mass and distance of the 
lens, which can potentially cause a detectable microlensing signature 
in direct imaging of Sgr A”. The values of mass and distance, which 
can potentially be probed with these possible baselines, are shown 
in fig. 11. The exact detectability of the microlensed shadow would 
require a detailed analysis of the visibility amplitude as a function 
of baselines which could distinguish the microlensed shadow and 
the true shadow with a given number of baselines, flux sensitivity, 
exposure time, and the dirty beam modeling of a VLBI facility. This 
will require a full numerical signal analysis for a given telescope 
specification which we leave for future exploration. This work will 
assume a simplified detection criterion of A > oq. 


4.1 Event duration 


Generically, we expect a relative lens motion relative to the true 
shadow. Therefore, the microlensed shadow changes its center, size, 
and shape with time. In the case of a moving lens, we then define the 
event duration as the duration for which the asymmetry is above the 
detection threshold og. For the relevant parameter space of lenses as 
depicted by fig. 11, we can estimate a sensitive region around Sgr A* 
within which a lens has to lie to create a measurable asymmetry 
by various possible baselines. We can then estimate the maximum 
threshold value of é by taking asymmetry in the large é limit, as 


R 
Tav? 


M 8.2kpcl-x 10-3 
26. ; 
6.5 mas 1; Mo D; x ok (36) 


&th 


Q 


OF 


Q 


Fig. 12 shows the variation in the duration of the microlensing 
event with the characteristics of the lens. The relative velocity of the 
lens is taken to be 100 km/sec with respect to the black hole shadow. 


Assuming the lens proper motion component along the radial di- 
rection of the true shadow is uy, the event duration of the microlens- 
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Figure 12. The above plot shows the duration of the microlensing event te 
during which asymmetry remains above 0.1% of R and its dependence on 
the lens mass M and the distance r from Sgr A* towards earth. The relative 
velocity of the lens is taken to be 100 km/sec with respect to the black hole 
shadow. 


ing event is given by, 


— 28th 
E= , 
Mr 
1/2 =3 
100 km/s { Ds M 10 
~20.6 yr (5 =] 1; “a — x) FJR" (37) 


It should be noted that to achieve the required contrast of inter- 
ference fringes in radio measurements, one needs to integrate the 
interference visibility amplitude over some minimum exposure time 
(depending upon the brightness of the source), and then repeated 
measurements over some time result in one black hole image. This 
collectively adds up to an effective imaging time tim, and hence, an 
observed black hole image will always be time-averaged over fim. 
Therefore, the duration for which the image of the shadow changes 
due to lensing should be larger than the imaging time. Otherwise, 
the microlensing feature would be washed out during the averaging. 
The imaging time used for Sgr A* in the current EHT observation is 
~ 1 day (Event Horizon Telescope Collaboration & Akiyama 2022a). 
For the potential space-based configurations, it has been estimated 
the imaging time of a month can achieve the microarcsec-level reso- 
lution (see ref. (Gurvits et al. 2022)), and in that case the microlensing 
events of durations above a month will be only observable. The exact 
imaging time will depend on many factors that go into implementing 
the space-based EHT-like facilities, which is beyond the scope of 
this paper to discuss in detail?. However, there is no fundamental 
limitation to reducing this imaging time, and it can be achieved with 
technological development. Hence, we will take a benchmark imag- 
ing time ft; = 1 day. We will thus assume that the target range of 
parameters of the lens, inducing lensing event with duration above 
one day (see e.g. dashed-line in fig. 12 drawn for the accuracy in ra- 
dial measurement to be 0.1% of R), will cause microlensing signals 


3 See ref. (Roelofs et al. 2019b; Gurvits 2020; Gurvits et al. 2021; Rud- 
nitskiy et al. 2022; Likhachev et al. 2022) for more details on the possible 
implementation of space interferometry. 
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which will not be washed out after averaging over the imaging time 
of a day. Furthermore, we can even expect a time-varying change 
in the shape of the shadow during the passes of a lens in front of 
Sgr A* with sufficient cadence and observational period of EHT-like 
facilities. 


4.2 Occurance of the microlensing due to stellar objects 


The shadow of Sgr A* is surrounded by various granular stellar 
objects. We estimate the probability that a microlensing event will 
be observed with a given baseline at any instant. This probability can 
be estimated to be 1 — e~* ~ t, where 7 is the optical depth given 
by (Kiraga & Paczynski 1994; Paczynski 1996), 


t= a ng PDD p 24D). (38) 
0 


The stellar density p( D7) towards the Galactic center of the Milky 
Way can be split into stellar bulge pp (r) and stellar disc pg(r). The 
stellar bulge can be modeled with the Hernquist (Hernquist 1990) 
density profile given by 


Mp a 
r) = ———, 39 
polr) on Ree (39) 
where r is the radial distance from the Galactic center, the mass 
parameter Mp = 2 x 1010 Mo (Valenti et al. 2016) and the scale 
radius a = 0.31 kpc (Li 2016). The stellar disc is described by 


Ma 


ela, (40) 
8nzor, 


Padisc(r) 


where My = 4.8 x 1010 Mo is the disc mass, and rg = 2.67 kpc is 

the disc scale length, and zg = 0.32 kpc is the scale height. Hence, 

the optical depth is found to be 
1073 


t œ 0.003—— | 


D 2 
sR - l (41) 


8.2 kpc 

We can also estimate the differential event rate of microlensing 
events using dI = dt/dte. In section 4.1, we established the sig- 
nificance of the duration te in a microlensing event, emphasizing its 
necessity to exceed tim. This ensures that the distinctive microlens- 
ing distortion signature remains discernible during the construc- 
tion of a single image using data collected within the time frame 
tim. Consequently, the total count of observable events within a 
specific baseline configuration relies on the distribution of events 
across different te values. To quantify this, we calculate the annual 
event count within twenty logarithmically spaced te bins spanning 
logjo(ti/1 yr) € [-2, 6]. The events in the i bins are counted using 


f n v p(Dı) 5 dD, 
AT; = f ia aes, (42) 
ti l dte 


where tf = ti X 100-4, 

In Figure 13, we present the distribution of microlensing event 
durations based on various accuracies in the asymmetry observable 
of Sgr A”. The prevailing trend indicates that a majority of events 
exhibit durations on the order of 107 days, 103 days, and 10+ days, 
corresponding to cg /R values of 1073, 1074, and 1075, respectively. 
Importantly, these durations far exceed the benchmark of one day 
(tim) and even a month for Earth-space configurations. As a result, 
the microlensing signatures of the majority of these events remain 
potentially observable and are not washed out during the imaging 
time. 
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Figure 13. This plot shows the distribution of the event durations te of 
the microlensing events potentially observable in the shadow of Sgr *. The 
velocity parameter is assumed to be 100 km/sec, and the mass M is 1 Mo. 
The histograms shown with red-dotted, green-dashed, and blue-sold lines 
correspond to the distributions assuming the uncertainties in the asymmetry 
observable relative to the radius of the true shadow (oa /R) to be 1073, 1074, 
and 1075 respectively. These are the typical accuracies achievable with the 
baseline configurations listed in the table. 1. The total number of events per 
year (T`) is also shown for each histogram. 


Hence, for a given baseline specified by the ratio ca /R, the total 
event rate can be estimated as 


r S” Zea 
- [5 


again assuming the abovementioned stellar density, 


3/2 = 
7 sine] Ds IMs 10-3 
T 14x 10™ y Tork (x3 =) y M FR: “A 


The subset of the above events which will induce a detectable magni- 
fication in size can then be obtained by dividing the above expression 


1/4 
by a constant factor of 38 (10-3/(0a/R)) and given by 


y D; \° 1Mo [1073 \!/ 
100 km/s \\ 8.2 kpc M \oq/R} ` 


(45) 


v PP) p 
iD; M 


Dj aD}, (43) 


L ~= 3.7x 107 yr! 


The above expression estimates the number of events potentially 
observed by baselines at the Moon and L2. These baselines in space 
will only have an event rate ~0.0014 event per year due to stellar 
granularity in the foreground of Sgr A*. Therefore, the stellar gran- 
ularity is insufficient to cause a microlensing effect on the shadow 
of Sgr A* very often. However, further investigation is needed with 
a numerical analysis of the visibility amplitude of the microlensed 
shadow, to properly model the accuracy of the radial measurements 
at various baselines. We plan to take this up in our future work. 


4.3 Potential enhancement in the event rate 


We observe infrequent microlensing events in the black hole shadow 
of Sgr A* caused by conventional stellar lenses toward the Galactic 
Center, even with sub-micro-arcsec resolution. In this context, we 
enumerate potential scenarios that could amplify the occurrence of 
microlensing events in Sgr A*. Additionally, we outline opportunities 
to detect microlensing distortions in the black hole shadows located 
at the cores of various galaxies. This encompasses not only the ev- 
ident case of M87, where the current EHT has already identified 
the shadow, but also numerous other shadows anticipated to become 
accessible through Earth-space interferometry with micro-arcsec or 
sub-micro-arcsec resolution in the future. 


4.3.1 Black hole cluster around the Galactic Center 


There has been a prediction of a cluster of about 20,000 stellar mass 
black holes within the central parsec of the Milky Way (Morris 
1993; Miralda-Escude & Gould 2000; Freitag et al. 2006). This 
cluster is formed as a result of dynamical friction, which predicts 
the density profile followed by these black holes will scale as rls 
and the average mass of these black holes is found to be about 
M =7 Mo (Miralda-Escude & Gould 2000). Although the cluster 
of black holes has not been detected yet, a density cusp of quiescent 
X-ray binaries, which are tracers of isolated black holes, has 
been detected (Hailey et al. 2018). This hints towards the existence 
of the cluster of black holes within a central parsec of the Milky Way. 


We provide a simple estimate of the contribution in the num- 
ber of microlensing events from this cluster of black holes, which 
can potentially be observable in the shadow of Sgr A*. In general, 
the density profile of the cluster of Nph number of black holes within 
ro radius from the galactic center can be written as, 


(46) 


5 NM [r \ 7 
p(r) z = ; 


167 i ro 


assuming each black hole has an average mass of M. The velocity 
dispersion of these black holes relative to the galactic center can be 
obtained using the Jeans equation as (Chaname et al. 2001) 


1 
68.5 km/s,/ —. (47) 


Following the formalism in the previous sec. 4.2, we cal- 
culate the event distribution over event durations under the 
assumptions Nph = 20,000 and rọ = 1 pc. The results are shown 
in fig. 14 under the assumption of various possible values of a/R 
achievable from Earth-space interferometry. We notice the presence 
of a cluster of black holes can provide a dominant contribution to the 
microlensing event rate in the shadow. In general, for Npp number of 
stellar mass black holes with average mass M within rọ radius from 
the galactic center, the total number of events can be summarized in 
the following expression 


fee Non (1pc\?( Ds \2 |7Mo 1073 
Dhe: ms 20,000 \ ro | \8.2 kpe M Ga/R’ 


(48) 


It is worth noting that the majority of these events exhibit 
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Figure 14. This plot shows the distribution of the event durations te of the 
microlensing events potentially observable in the shadow of Sgr A* assuming 
the black hole mass M to be 7 Mo. The total number of events per year (T) 
is shown for each histogram. However, the event rate of durations above | day 
are 7 x 1074, 5 x 1073, and 4 x 1072 corresponding to Ca/R to be 1073, 
1074, and 1075 respectively. These are the typical accuracies achievable with 
the baseline configurations listed in table. 1. 


durations below the benchmark imaging time fim = 1 day. 
Specifically, event rates with durations exceeding one day amount 
to 7 x 1074, 5x 1073, and 4 x 1072, corresponding to a/R 
values of 1073, 1074, and 1075 respectively. Thus, detecting these 
short-duration microlensing events remains challenging due to 
the potential clustering of black holes around Sgr A*. However, 
the considerable number of events with durations smaller than 
tim might induce a blurring effect during shadow imaging. This 
microlensing-induced blurring could potentially be detected by 
reaching an irreducible angular scale of the order of the Einstein 
angle. Below this scale, the image may remain unresolved even with 
the capability of a VLBI facility. 


4.3.2 Other shadows 


We have illustrated the impact of microlensing in the shadow of 
SgrA* and examined the potential observability of this phenomenon 
under various circumstances. A comparable investigation could be 
conducted for the shadow of M87, evaluating the feasibility of ob- 
serving such an event by considering a well-informed distribution 
of compact objects toward M87. Additionally, as black hole imag- 
ing achieves higher resolution, numerous other shadows within the 
capability of EHT-like facilities will be detectable. We anticipate 
that sub-millimeter interferometric observations with approximately 
0.1yas resolution and ~ 1 Jy sensitivity may enable the observa- 
tion of over ~ 10° supermassive black hole shadows (Pesce et al. 
2021). Consequently, even if the event rates for individual shadows 
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are low, the increased number of measured shadows can enhance the 
likelihood of observing a microlensed shadow. 


5 SUMMARY 


We have presented a formulation of the phenomenology of microlens- 
ing via observing directly imaged shadows which are potentially de- 
tectable in a near-future EHT-like VLBI placed in space. We found 
an exact analytic expression of the shift in the center and the radial 
profile of the microlensed shadow. The shift in the center, magnifica- 
tion in the average size, and the asymmetry in shape (defined as the 
standard deviation in the radial measurements along various position 
angles) are the key signals imprinted on a microlensed shadow. The 
maximum average radius of the microlensed shadow could reach 
up to 3R/2, and the corresponding maximum asymmetry would be 
R/8V2. 

In the case of the microlensed shadow of Sgr A*, we find that the 
maximum average radius is 36.535 jas, and the maximum asymme- 
try can reach up to 2.152 uas. We presented a range of the parameter 
space of lenses that can be probed with the phenomena of microlens- 
ing via future baselines in space. The event rate for future baselines 
at the Moon and Ly is found to be 1.4 x 1073 yr! for lenses of mass 
~ Mo in the stellar component. Therefore, observing the microlens- 
ing effect due to such lenses will always be a rare effect. However 
the signal would be boosted by an anticipated cluster of stellar mass 
black holes around the galactic center that would enable us to detect 
the effect of microlensing in the shadow of Sgr A*. 

It is worth mentioning here that a more exotic but currently vi- 
able scenario of primordial black holes (PBHs) as dark matter (Carr 
& Kuhnel 2021) with exciting prospects for future detection via 
gravitational wave signatures of solar or sub-solar mass candidates 
(Bagui et al. 2023), should imply an even higher black hole density 
towards the Galactic Center. This is anticipated in the possible pres- 
ence of a dark matter spike within the gravitational influence radius 
(~ 0.2 pc) of Sgr A*. This hypothesis is motivated by theory (Gon- 
dolo & Silk 1999), constrained by stellar orbits (Shen et al. 2024), 
and, less directly, by the observed old stellar density distribution 
(Habibi et al. 2019). Hence the PBH interpretation of dark matter 
may eventually be tested by estimating the microlensing event rate 
due to PBHs and comparing it with future high-resolution imaging of 
the shadow. Additionally, microlensing black hole shadows with the 
future space-based VLBI facilities will open up exciting possibilities 
for observing other shadows and hence probing the existence and 
properties of compact objects in the mass range of 1076 — 104 Mo 
towards the galactic centers. 

Lastly, even though we have demonstrated the effect of microlens- 
ing on the boundary of the shadow, the formalism presented in this 
work is equally valid for an impact on the predicted gravitational 
lensing rings. Microlensing will make the lensing rings larger, asym- 
metric, and shifted as compared to their true shapes. 
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APPENDIX A: OTHER OBSERVABLE PLOTS 


We show the induced shift in the lensed center for different angular 
separations of the lens from the shadow of Sgr A* (fig. Al) as a 
function of the lens’s mass and distance from the center of Sgr A*. 

The induced magnification is derived in fig. A2 for the average 
radius of the microlensed shadow and different angular separations 
of the lens from the center of Sgr A*. 

Finally, we depict the induced asymmetry in the microlensed 
shadow corresponding to different angular separations of lenses of 
different masses from the center of Sgr A* (fig. A3). 
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Figure A1. This figure depicts the lens’s mass and distance from the center of Sgr A* and the corresponding induced shift in the lensed-center. Various panels 
correspond to the different angular separation of the lens from the shadow of Sgr A*. 
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Figure A2. This figure depicts the lens’s mass and distance from the center of Sgr A* and the corresponding induced magnification in the average radius of the 
microlensed shadow. Various panels correspond to the different angular separation of the lens from the shadow of Sgr A*. 
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Figure A3. This figure depicts the lens’s mass and distance from the center of Sgr A* and the corresponding induced asymmetry in the microlensed shadow. 
Various panels correspond to the different angular separation of the lens from the shadow of Sgr A*. 
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